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1 Radiation-hydrodynamics modeling — overview 

In this paper we discuss three issues in the context of three-dimensional (3D) 
hydrodynamical model atmospheres for late- type stars, related to spectral 
line shifts, radiative transfer in metal-poor 3D models, and the solar oxygen 
abundance. To establish the context we start by giving a brief overview about 
the model construction, taking the radiation-hydrodynamics code C0 5 BOLD 
(Conservative COde for the Computation of COmpressible COnvection in 
a BOx of L Dimensions with L=2,3; [3]) and the related spectral synthesis 
package Linfor3D as examples. 

Based on a Godunov-type finite volume approach, C0 5 BOLD provides 
the time-dependent solution for a one-component compressible radiating fluid 
in an external gravity field on a fixed, non-staggered 3D Cartesian grid (allow- 
ing variable spacing). Operator splitting separates Eulerian hydrodynamics, 
optional tensor viscosity, and radiation transport. Directional splitting de- 
composes the 3D hydrodynamics problem into ID sub-steps which are treated 
by an approximate Riemann solver of Roe type, modified to work with an 
arbitrary equation of state, and to properly handle an external gravity field. 
This scheme is very robust and well adapted to handle transonic flows and 
shocks in a highly stratified medium. By design, the code guarantees the 
numerical conservation of mass, momentum, and energy. For any prescribed 
chemical composition, C0 5 BOLD uses a tabulated equation of state taking 
into account partial ionization of HI, He I, and Hell, as well as the formation 
and dissociation of H2 molecules. 

The role of radiation in the hydrodynamical simulations is to describe the 
energy balance due to radiative heating and cooling. The radiative energy 
exchange is computed from the solution of the non-local transfer equation on 
a system of a large number of rays traversing the computational volume un- 
der different azimuthal and polar angles. Realistic stellar opacities are used, 
optionally based on ATLAS or MARCS opacity data. The frequency depen- 
dence of the radiation field is treated in a multi-group approximation - the 
so-called opacity binning method (OBM; [9, 7, 10]) - where frequencies are 
sorted into a small number of bins (typically 4 ... 6) according to the ratio 
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Fig. 1. Left panel: Time series of line profiles of a Fel line at 6082 A in a 3D 
solar CO'BOLD model. The wavelength is given as corresponding Doppler velocity 
with respect to the line's laboratory wavelength. The line profiles have been offset 
proportional to time running from top to bottom. The time interval between profiles 
is the same, fluctuations in the continuum brightness cause the non-equidistant 
appearance. Right panel: Contour plot of correlation coefficients of the intensity 
between two wavelength points in the profiles shown in the left panel. 

of monochromatic to Rosseland optical depth. So far, strict LTE is assumed, 
thus scattering cannot be treated. Radiation pressure is ignored. 

The code Linfor3D accepts C0 5 BOLD models as background structures 
on which spectral synthesis calculations at high wavelength resolution - usu- 
ally focusing on one particular spectral line - can be performed. When calcu- 
lating the emergent spectrum, Linfor3D takes into consideration the full 3D 
flow geometry including Doppler shifts caused by the macroscopic hydrody- 
namical velocities. It represents the effects of thermal and pressure broaden- 
ing in standard fashion, but leaves out the ad-hoc broadening mechanisms of 
micro- and macro-turbulence introduced in ID atmospheric models. Similar 
to C0 5 BOLD, strict LTE is assumed in Linfor3D. Resulting spectral line pro- 
files provide detailed information about intrinsic line shapes, and convective 
line shifts with respect to a line's laboratory wavelength. 

2 High precision line shifts from 3D models? 

A C0 5 BOLD simulation constitutes a statistical realization of the atmo- 
spheric flow field in the stellar surface layers. If one is not interested in 
studying time-variable phenomena but only in the mean state of the atmo- 
sphere, the fluctuations present a noise source which limits the precision to 
which flow and related spectroscopic properties can be determined. This is 
similar to the observational situation where the intrinsic variability of a star 
limits the precision to which its radial velocity, e.g. in planet searches, can 
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be measured. The left panel of Fig. 1 illustrates the temporal variability of a 
Fe I line calculated from a hydrodynamical solar model. Each line profile is a 
horizontal average over the surface of the computational box. Shown are 25 
instants in time which are sufficiently separated that they can be considered 
statistically uncorrelated. We ask: what is the precision due to the statistics 
(ignoring systematic effects) to which we can determine the line shift? 

It is straight forward to show (see [8]) that the expectation value of the 
disk-integrated line profile corresponds to the expectation value of the profile 
of the local hydrodynamical model. Hence, the statistical uncertainties of the 
profile emerging from the model directly correspond to the uncertainties of 
the predicted disk-integrated profile. From Fig. I it is obvious that the the 
statistical fluctuations are not just pixel-to-pixel random noise like in the case 
of photometric Poisson noise. The line profiles change their overall shape, i.e. 
different wavelength points show a considerable degree of correlation. The 
linear correlation coefficient between intensities at two wavelength points 1 
and 2 is given by 



I denotes the intensity, (.} the temporal average. Ali = /» — {h) is the inten- 
sity deviation from the mean. The right panel of Fig. I shows the correlation 
matrix of the example line depicted. 

In order to quantify the line shift A s we need a model of the procedure by 
which it is measured, which in turn emerges from the definition of A s . Here, 
we assume that the measuring procedure of A s can be described by a function 
A of potentially all available (assumed discrete) intensities If. A s = A (7j). In 
order to make algebraic headway we simplify and linearize A around the 
expectation value of the line profile described by the values To leading 
order in AI we obtain for the variance of the line shift the standard expression 
of the error propagation for correlated variables 



The summation is performed over all pixels which are relevant for the mea- 
surement of A s . Equation (2) emphasizes the role of the covariance matrix 
of the intensities {Ali Alj ) - or equivalently the standard deviations of the 
intensities and their correlation matrix - plays for the magnitude of the un- 
certainty of the line shift. In the present context we discussed line shifts but 
relation (2) of course also holds for other measures like, e.g., the equivalent 
width of a line. The statistical quantities in relation (2) can be estimated 
from the time series provided by the hydrodynamical model. Asymptotically, 
for a given line one will arrive at a fixed value for the correlation matrix 
C[Ii,Ij]. If one wants to improve the accuracy of the line shift one has to 



C[h,h] = 



{Alj Alj) _ {hh) - {h){h) 



(1) 
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beat down the uncertainties in the intensities ai i . This can be achieved by 
longer time series or larger horizontal extent of the hydrodynamical model. 

Our example Fel line shows a RMS temporal scatter of its position of 
0.16kms _1 . The value was obtained by directly (and somewhat heuristi- 
cally) measuring the location of the line core without formalizing the process 
by explicitcly constructing a measurement function A. The statistical in- 
dependence of the 25 individual snapshots implies an uncertainty of about 
30 ms" 1 for the line shift. While the specific value depends on the chosen 
line and selection of snapshots we think that it gives an indication of the 
precision one is typically working with in todays hydrodynamical standard 
models. Higher precision is possible but computationally also more costly. 
Of course, at some point real uncertainties will be dominated by systematic 
shortcomings of a model. 

3 3D radiative transfer in metal-poor atmospheres 

As mentioned earlier, the radiative transfer in the 3D models is commonly 
approximated by the opacity binning method (OBM) assuming strict LTE. 
While the approach is working fine in atmospheres of about solar metallicity, 
metal-poor atmospheres pose a challenge to the OBM. At first glance, this 
may come as surprise because the dramatic decrease of the number of spectral 
lines relevant for the radiative energy exchange should simplify the radiative 
transfer. However, the actual situation is quite different. First, scattering in 
the continuum becomes important for the thermal structure of metal-poor 
atmospheres. In the OBM, scattering is treated as true absorption so that 
one must expect some effects on the resulting temperature structure. Second, 
experience has shown that the OBM does not work as accurately in metal- 
poor atmospheres as in atmospheres of solar metallicity It turned out that 
this deficit is not related to the treatment of the line blocking but already 
shows up for the radiative transport in the continuum. 

Figure 2 shows an example of an atmosphere of a metal-poor giant. Plot- 
ted are temperature profiles of ID ATLAS6 (see [6]) model atmospheres in 
radiative-convective equilibrium. The only difference among the models is 
the way in which the radiative transfer was treated. In three cases labeled 
"scattering", "no scattering", "scattering as true absorption" a high wave- 
length resolution was employed, and scattering was treated exactly, scatter- 
ing opacity was neglected, or treated as true absorption, respectively. The by 
far dominating scattering opacity under the considered conditions is Raylcigh 
scattering by hydrogen atoms. As evident from Fig. 2, the temperature struc- 
ture is noticeably influenced by scattering. The OBM used in the 3D models 
was also implemented in the ID atmosphere code and a resulting radiative- 
convective equilibrium calculated. Comparison with the exact radiative trans- 
fer solution shows a close correspondence from the deep layers up to lower 
optical depth of logr « —3. However, while useful in practice this is only 
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Fig. 2. Comparison of ID model atmospheres (T e s = 5000 K, log(g) = 2.94, 
[M/H] — —2) in radiative-convective equilibrium based on different treatments 
of the radiative transfer, (for details see text) 

fortuitous. The OBM based temperature structure should in fact follow the 
"scattering as true absorption" case since in the OBM scattering opacity is 
treated as true absorption. At present, the reason for the rather poor per- 
formance of the OBM is unclear. Identifying its cause, improving the OBM, 
and including scattering effects in 3D geometry arc challenges to be met in 
near-future developments of 3D model atmospheres for metal-poor stars. 

4 3D models and the solar oxygen abundance 

Recent spectroscopic abundance determinations by [1], based on a 3D hy- 
drodynamic model atmosphere, led to a much debated downward revision of 
the solar C, N, and O abundances. Their result for the oxygen abundance is 
logeo =8.66 ±0.05 (on the scale log en = 12), causing a dramatic deteriora- 
tion of the agreement between the thermal structure derived from helioseis- 
mic inversions and theoretical solar models, respectively. Motivated by this 
problem, we (see [2]) are currently using a 3D C0 5 BOLD simulation with 
5-bin frequency-dependent radiative transfer based on MARCS opacities to 
see whether the results by [1] can be confirmed. This independent redetermi- 
nation of the solar oxygen abundance is based on 2 forbidden and 7 permitted 
OI lines, using a number of different observations, including both disk-center 
("intensity") and full-disk ("flux") spectra. In addition to 25 snapshots from 
the simulation, we also derive abundances from different ID atmospheres for 
comparison. Special care is taken to provide realistic error estimates. 

The following preliminary conclusions can be drawn at this point: (i) 
"intensity" and "flux" spectra give practically the same result, (ii) the oxygen 
abundance derived from the 3D C0 5 BOLD simulation is only slightly lower 
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(by -0.04 dex) than that derived from the ID empirical model by [4] (hereafter 
HM), indicating that the 3D mean model and the ID HM model have very 
similar temperature structures in the relevant layers, (iii) the unknown cross 
sections for neutral particle collisions introduce uncertainties in the NLTE 
corrections for the O I triplet lines of up to 0.1 dex; depending on the weight 
of these lines, the resulting error in the mean oxygen abundance derived 
from our set of lines is about 0.05 dex. (iv) Our preliminary best estimate 
for the solar oxygen abundance is logeo = 8.72 ± 0.06, which is close to the 
value recommended by [5], logeo = 8.736 ± 0.078. A remaining problem of 
our analysis is that the two forbidden [OI] lines give significantly different 
abundances, which cannot be explained by NLTE-cffccts or deficiencies of 
the model atmosphere. We hope to resolve this problem by analyzing the 
observed center-to-limb variations of these two line profiles. 

5 Remarks on precision spectroscopy and 3D models 

Hydrodynamical model atmospheres arc on their way of becoming a stan- 
dard tool for the analysis of stellar spectra. Their ability of making detailed 
predictions about the shape of spectral lines in convective atmospheres can 
only fully exploited if observed spectra of sufficient resolution are available. 
Ideally, spectrographs should be able to provide a spectral resolutions above 
10 5 - something that we would like instrument builders to keep in mind. 

High-fidelity abundance work benefits from the theoretical knowledge of 
the precise intrinsic line shape. However, in practice one is nonetheless often 
confronted with ambiguities, e.g. in the case of blends, which remain unre- 
solved by considering disk-integrated line profiles only. The center-to-limb 
variation of a line shape can provide crucial further constraints. Combining 
interferometry with high-resolution spectroscopy (like in the UVES-I project 
presented by A. Quirrenbach, this volume) can open-up this source of infor- 
mation for stellar work. 
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